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We study the impact of a stochastic background of primordial magnetic fields on the scalar contri- 
bution of CMB anisotropics and on the matter power spectrum. We give the correct initial conditions 
for cosmological perturbations and the exact expressions for the energy density and Lorentz force 
associated to the stochastic background of primordial magnetic fields, given a power-law for their 
spectra cut at a damping scale. The dependence of the CMB temperature and polarization spectra 
on the relevant parameters of the primordial magnetic fields is illustrated. 
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I. INTRODUCTION 

Large scale magnetic fields are almost everywhere in 
the universe, from galaxies up to those present in galaxy 
clusters and in the intercluster medium [l| . The origin of 
these magnetic fields depends on the size of the objects 
and may become mysterious for the largest ones. The 
dynamo mechanism provides a mechanism to explain the 
observed magnetic field associated to galaxies, whereas 
those associated to clusters may be generated by gravi- 
tational compression starting from an initial seed. 

The requirement of an initial seed for magnetic fields 
observed in galaxies and galaxy clusters leads directly to 
question the existence of primordial magnetic fields in the 
early universe. Cosmology described by an homogeneous 
and isotropic expanding metric neither supports a uni- 
form magnetic field nor a gravitational amplification of 
gauge fields because of conformal invariance; the genera- 
tion of large scale magnetic fields has therefore generated 
a lot of interest. A stochastic background (SB) of primor- 
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dial magnetic fields (PMF) can provide the initial seeds 
for the large-scale magnetic fields observed and can leave 
imprints on different observables, as the CMBpattern of 
temperature and polarization anisotropies 0, y[ and the 
matter power spectrum. 

A SB of PMF carries zero energy and pressure at ho- 
mogeneous level in a Robertson- Walker metric. It carries 
however perturbations, of any kind, i.e. scalar, vector 
and tensor, and it is usually studied in a quasi-linear ap- 
proximation, i.e. its EMT - quadratic in the magnetic 
field amplitude - is considered at the same footing as 
first order terms in aperturbative series expansion. Vec- 
tor and tensor 0)0)0 me tric perturbations sourced 
by a PMF SB have been object of several investigations; 
beyond the technical simplicity of vector and tensor over 
scalar, a perfect fluid cannot support vector and tensor 
perturbations at linear order and therefore represent a 
key prediction of a PMF SB. We know however that tem- 
perature and polarization anisotropies sourced by scalar 
fluctuations with adiabatic initial conditions are a good 
fit to the whole set of observations; it is therefore crucial 
to investigate how a PMF SB can modify these scalar 
fluctuations. Analytic [1] and numerical [1, QiiGjl works 
in this direction have already been made. However a 
detailed analysis which takes into account the Lorentz 
force on baryons, a careful treatment of initial conditions 



2 



and an accurate treatment of the Fourier spectra of PMF 
energy-momentum tensor is still lacking. As is clear in 
the following, our work address carefully both these is- 
sues. 

The goal of this paper is to investigate the impact 
of a stochastic background (SB) of primordial magnetic 
fields (PMF) on scalar cosmological perturbations and in 
particular on CMB temperature anisotropics and matter 
power spectrum. Our paper is organized as follows. In 
Section II we review how to add a fully inhomogenous SB 
of PMF treated in the one-fluid plasma description [2j to 
the Einstein-Boltzmann system of equations. In Sections 
III and IV we review the baryons evolution and we give 
the initial conditions for cosmological perturbations in 
a form suitable to be plugged in most of the Einstein- 
Botzmann codes. In Section V we give the PMF energy 
density and Lorentz force power spectra and compare our 
results with the ones given in the literature. In Sections 
VI- VIII we show the results obtained by our modification 
of the Einstein-Boltzmann code C AMB [l8[ for cosmolog- 
ical scalar perturbations, CMB spectrum of temperature 
and polarization, matter power spectrum, respectively. 
In the Appendix we show the detailed calculations for 
the convolution integrals leading to the energy density 
and Lorentz force, starting from a power-law spectrum 
sharply cut at a given scale for the PMF. 



II. STOCHASTIC MAGNETIC FIELDS AND 
COSMOLOGICAL SCALAR PERTURBATIONS 

We model a SB of PMFs as a fully inhomogenous com- 
ponent, considering B 2 at the same level of metric and 
density fluctuations in a perturbative expansion *. Al- 
though a SB of PMFs carries no energy at the homoge- 
neous level, it affects scalar cosmological perturbations in 
three different ways. First, inhomogeneous PMFs carry 
energy density and pressure and therefore gravitate at 
the level of perturbations. Second, inhomogeneous PMFs 
have anisotropic stress - differently from perfect fluids - 
which adds to the photon and neutrino ones, with the 
caveat that the photon anisotropic stress is negligible be- 
fore the decoupling epoch. Last, but not least, the in- 
duced Lorentz force acting on baryons, affects also pho- 
tons during the tight coupling regime. 

Since the EMT of PMF at homogeneous level is zero, 
at linear order PMFs evolve like a stiff source and there- 
fore it is possible to discard all the back reactions of the 
fluid or gravity onto the SB of PMF. Before the decou- 
pling epoch the electric conductivity of the primordial 



Note that in such a way we do not take into account the modifica- 
tion of the sound speed of baryons induced by PMFs, pionereed 
in |12| , since it would be technically of second order in the equa- 
tions of motion. However, since the baryons speed of sound goes 
rapidly to zero in the matter dominated era, this effect, leading 
to a shift in the Doppler peaks, may be anyway important. 



plasma is very large, therefore it is possible at the first 
order to consider the infinite conductivity limit. In this 
limit the induced electric field is zero. Within the infinite 
conductivity limit the SB of PMF time evolution simply 
reduces to : B(x, t) = B(x)/a(r) 2 . t 

The evolution of the metric perturbations in the pres- 
ence of PMF is governed by the Einstein equations: 

G^=87r(T^+T™ F ), (1) 

In the approximation in which the induced electric field is 
vanishing (i.e. the infinite conductivity limit) the energy 
momentum tensor of the electromagnetic field becomes: 
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In the Fourier spaced the Einstein equations with the 
contribution of PMF in the synchronous gauge are: 

k 2, q - ^Hh = 47rGa 2 (E„ p n S n + p B ) ■ 

k 2 f] = 4 7 rGa 2 E„( / 9„ +P„)0„, 
h + 2Hh-2k 2 r] = -8irGa 2 (Z n c 2 snPn 5 n 



h + 6ij + 2H(h + 6i))-2k 2 r) = -24ttGo 2 



where by n we mean the number of components, i.e. 
baryons, cold dark matter (CDM), photons and neutri- 
nos. The conservation of the PMF EMT - V /1 r^ PMF = 
- simply reduces to : 



Pb , T 
v B = — +L. 



(7) 



where a B represents the PMFs anisotropic stress and L 
the Lorentz force. The energy density of PMF evolves 
like radiation: /?b(x, r) = ps(x, ro)/a(r) 4 . 



III. BARYONS EVOLUTION 

The presence of PMFs in a plasma which contains 
charged particles induces a Lorentz force on these par- 
ticles, that, in the primordial plasma, are baryons. The 



t We choose the standard convention in which at present time to, 
a(t ) = 1. 

t As Fourier transform and its inverse, we use - in agreement with 
E1-: 



(2tt)3 

where Y is a generic function 



Y(x,r), Y(x,t) = J d ke Y(k,r) . 

(5) 
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general expression for the Lorentz force is 



IV. INITIAL CONDITIONS 



Li(x,T ) 



1 

4tt 



1. 



B j (x)V i B i (x) - -V t £H(x) 



(8) 



where L(a;,r) = L (^ T0 ) . 

We are interested only in the scalar perturbations and the 
scalar part of the Lorentz force defined as V 2 L^ = VtLi 
is therefore: 



V 2 L^ = — 
An 



(V 4 B J (x))V J S l (x)-iv 2 i? 2 (x) . (!)) 



In the presence of an electromagnetic source the con- 
servation equations of the baryon component of the pri- 
mordial fluid becomes: 



V AT^ baryons oc F^J U 



(10) 



where is the quadrivector of the density current and 
is the Maxwell tensor. The primordial plasma can 
be considered globally neutral, this leads to Jo = 
and therefore to the fact that the energy conservation 
of baryons is not modified by the presence of the Lorentz 
term. The Euler equation for baryons in instead affected 
by the Lorentz force and the scalar part is therefore 0]: 



= -W b + k 2 c 2 b 5 b - k 
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Now we study how the tight-coupling regime is modified 
by the presence of a SB of PMF [14| . The Euler equation 
for photons during the tight-coupling regime is: 



= k 



+ an e OT\pb 
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Combining the photons and baryons equations gives: 



9 b = 
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with: 

6b — 6 7 




The photon Euler equation in tight coupling regime 
instead is: 



-R-Htib + Wb- c 2 k 2 S b + k 2 — 



(13) 



We note that there is a term depending on the Lorentz 
force which disappears when the tight coupling ends, 
leaving the normal Euler equation for the photon veloc- 
ity. 



In order to study the effect of a PMF SB on scalar 
cosmological perturbations, the initial conditions for the 
latter deep in the radiation era are required (see for the 
results in the longitudinal gauge). The magnetized adi- 
abatic mode initial conditions in the synchronous gauge 
are given by [151 ] : 
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where R v — p v /{p v + p 7 ) and C\ is the constant which 
characterize the regular growing adiabatic mode as given 
in [l3j]. We have checked that the result reported in [ll| 
and ours [l5[ agree. 

Note how the presence of a SB of PMFs induces a new 
independent mode in matter and metric perturbations, 
i.e. the fully magnetic mode. This new independent 
mode is the particular solution of the inhomogeneous 
system of the Einstein-Botzmann differential equations: 
the SB of PMF treated as a stiff source acts indeed as 
a force term in the system of linear differential equa- 
tions. Whereas the sum of the fully magnetic mode with 
the curvature one can be with any correlation as for an 
isocurvature mode, the nature of the fully magnetic mode 
- and therefore its effect - is different: the isocurvature 
modes are solutions of the homogeneous system (in which 
all the species have both background and perturbations), 
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whereas the fully magnetic one is the solution of the in- 
homogeneous system sourced by a fully inhomogeneous 
component. 

It is interesting to note the magnetic contribution 
drops from the metric perturbation at leading order, al- 
though is actually larger than the adiabatic solution for 
photons, neutrinos and baryons (the latter being tightly 
coupled to photons deep in the radiation era). This is 
due to a compensation which nullifies the sum of the 
leading contributions (in the long-wavelength expansion) 
in the single species energy densities and therefore in the 
metric perturbations. A similar compensation exists for 
a network of topological defects, which does not carry 
a background energy-momentum tensor as the PMF SB 
studied here §. 



where ey/ is the totally antisymmetric tensor, Pb and 
Ph are the non-helical and helical part of the spectrum 
for the B amplitude, respectively. Scalar cosmological 
perturbations only couple to the non-helical part of the 
spectrum and we shall therefore consider only Pb in the 
following. 



A. Magnetic Energy Density 



V. MAGNETIC FIELD POWER SPECTRA 

Power spectra for the amplitude and the EMT of 
SB of PMF have been subject of several investigation 
0, S @, [§] ■ We shall work in the Fourier space accord- 
ing to Eq. ©• We shall consider PMFs with a power 
law power spectrum, which therefore are characterize by 
two parameters: an amplitude A and a spectral index 
ns ■ PMFs are suppressed by radiation viscosity on small 
scales: we approximate this damping by introducing an 
ultraviolet cut-off in the power spectrum at the (damp- 
ing) scale kn- 

The two-point correlation function for a statistically 
homogeneous and isotropic field is 



(B?(fc)^(A : '))=5 3 (fc-fc') 



(Sij kikj) - h 



tiil\p H {k) , (15) 

Note however that a network of topological defects jiocs not scale 
with radiation and interacts only gravitationally with the rest of 



As is clear from Eqs. (1-3), the EMT for PMF is 
quadratic in the field amplitude. The PMF energy den- 
sity spectrum is 0: 



\PB{kt 
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(16) 

As for the two- 



p\k-p\ y / {k' 2 +p 2 -2kpcos9) 

point function in the coincidence limit, for several physi- 
cal spectral indexes and PMF configurations such convo- 
lution is not finite. There are in general problems both 
on large and short scales. Since the spectrum of the com- 
ponents of PMF EMT are relevant for the final impact 
on cosmological perturbations and CMB anisotropics, it 
is better to address this point in much more detail with 
respect to what is present in literature. 

matter, i.e. a Lorcntz term is absent. 
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The usual choice in the literature is to modify the scalar part two-point function of Eq. (|15p for zero helicity as 

/ <S 3 (k - k' ) (Sa - k kj ) f or k < k D 



(B*(k)B 3 (k')) 



\0 



for k > kj 



r 



with 



P B (k) = A[- 



(17) 



where fc* is a reference scale. With such choice the two- 
point function in the coincident limit (the mean square 
of the magnetic field) is: 
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(18) 



It is also usual in the literature to give the amplitude of 
B at a given smearing scale k$ by imposing a Gaussian 
filter : 



(B 2 (x)) ks = d 3 kP B (k)e' k2 / k s 
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FIG. 1: Plots of the different ways of computing the magnetic 
power spectra (in units of 4nAk D B+3 /k" B ) versus ub for ks = 
k D /2 (left) andfcs = k D (right). {B 2 }, (B 2 ) ks , (B 2 ) cut , ks are 
represented by solid, dotted and dashed lines, respectively. 



diveregence in k in the exact result; for n B < —3/2 the 
spectrum is no more white noise for k << fcjj. Fig. @ 
shows the dependence of k 3 \p B (k)\ 2 on ng at fixed (B 2 ). 

Our result are different from the one reported in the 
literature ISl, which is 
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By smearing the magnetic power spectrum and integrat- 
ing for k < kn, one gets: 



d 3 kP B {k)e- k2/k " 
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where the incomplete Gamma function T(... , ...) [l7j has 
been introduced. Note how n B > —3 in order to prevent 
infrared divergencies either in the mean square field or 
the amplitude of the field smeared at a given scale. In 
the following by (B 2 ) we mean the value given by Eq. 
(fT8|) . Fig. shows how (B 2 {x))k s may be much larger 
than (B 2 ) for n B > 0. 

The exact result for the Fourier convolution leading to 
the magnetic energy density Fourier square amplitude is 
one of the new main results of this paper. The convolu- 
tion involves a double integral, one in the angle between 
k and p and one in the modulus of p. The integral in 
the angle, often omitted in the literature, is the reason 
for having the result for |ps(fc)| 2 non vanishing only for 
k < 2ko. The detailed calculations for the energy density 
convolutions are given in Appendix A for several values of 
n B . The generic behaviour for k << k B and n B > —3/2 
is white noise with amplitude 



\PB(k)f 



16irk 2n (3 + 2n B ) 
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and then goes to zero for k — 2ko, which is a result 
obtained by performing correctly the integral. The pole 
for n B = —3/2 in Eq. (f2"Tj) is replaced by a logarithmic 
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FIG. 2: Plot of magnetic energy density power spectrum 
k 3 \p B {k)\ 2 in units of (B 2 ) 2 /(1024tt 3 ) versus k/k D for dif- 
ferent ub for fixed (B 2 ). The different lines are for ub = 
—3/2,-1,0,1,2,3,4 ranging from the solid to the longest 
dashed. 




FIG. 3: Comparison of magnetic energy density convolution 
k%\pB{k)\ 2 obtained in this paper (dotted, solid) in units 
of (B 2 ) 2 /(1024tt 3 ) and the one in Eq. ([22} (dashed, long- 
dashed) versus k/k D for n B = 2, -3/2 with fixed (B 2 ). 



and is not limited in k. In Fig. §5§ we show the difference 
between the literature result [f| and our result for n B = 
2,-3/2. 



B. Lorentz Force 



As is clear from previous sections, we also need the Lorentz force 



|L(fc)|2 = 128^8 / d3p Pb ^ Pb ^ - P'^ 1 + ^ + - V® ' 



(23) 
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and the magnetic anisotropic stress 
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where 7 = fc • p, = k ■ (fe — p)/(k\k — p\) and = p ■ (k — p )/(p\k — p\). 





FIG. 4: Plot of the Lorentz force power spectrum k 3 \L(k)\ 2 
in units of (_B 2 ) 2 /(10247r 3 ) versus k/ko for different ns for 
fixed (B 2 ). The different lines are for iib = —3/2, —1, 0, 1, 2, 3 
ranging from the solid to the longest dashed. 



FIG. 5: Comparison of the magnetic energy density and 
Lorentz force power spectra versus k/ko for fixed (B 2 ). The 
solid (medium dashed) and long-dashed (short-dashed) lines 



are respectively for 
(n B = -3/2). 



fc^p s (A:)| 2 and k J \L(k)\ 2 for n B 



We decide to compute the spectrum of the Lorentz 
force and obtain the anisotropic stress by Eq. (|7|). The 
exact computation for the Lorentz force power spectrum 
is given in Appendix B for several values of tib. A 
term —ps can be easily identified in Eq. (0); since we 
know from the exact computation that the integral of 
Pb(p) Pb{\^ — p I ) (1 + /i 2 ) is larger than the remaining 
piece in Eq. (j23|) we chose the the signs for ps(fc) and 
L(k) as opposite. 

Fig. ((4| shows the dependence of A: 3 |L(fc)| 2 on at 
fixed (P> 2 ). Fig. ([5]) compares the approximation L(k) ~ 
— ps(fc) suggested in Ref. [1] with the exact calculation. 
As can be checked in Appendix B, our exact calculations 
for the values of ub studied here show that 
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\pB{k)\ 2 for k « k D . 



(25) 



effect on the baryon velocity is crucial. Much later than 
the larger between the decoupling time and the time at 
which sound speed of baryons is effectively zero, baryons 
velocity can be approximated as: 



nlatc 



La\ t 
Pb ) a 



(26) 



during the matter dominated era. Our modified Einstein- 
Boltzmann code reproduces correctly this asymptotic 
regime for different wavelengths, as can be seen by Fig. 
([6]). The corresponding effects on the density contrasts 
for the same wavelengths are shown in Fig. {?]). In Fig. 
© the effects due to the pure magnetic mode and due 
to the correlation with the adiabatic mode are shown. 
Fig. ^ displays the importance of the Lorentz term 
compared to the purely gravitational effect. 



VI. RESULTS FOR COSMOLOGICAL 
PERTURBATIONS 

In order to study the effects of a SB of PMFs on CMB 
anisotropics and matter power spectrum we modified the 
CAMB Einstein-Boltzmann code [18| (June 2006 version) 
by introducing the PMF contribution in the Einstein 
equations, in the evolution equation for baryons and ini- 
tial conditions, along Eqs. (|6I26I14[) . 

We note that implementing the baryons evolution as 
from Eq. (|26|) , the MHD approximation in a globally neu- 
tral plasma is used up to the present time: this makes the 
Lorentz term non-vanishing up to the present time. Al- 
though the term L/ p\, in Eq. (|2"6")) decreases with time, its 



VII. RESULTS FOR CMB TEMPERATURE 
AND POLARIZATION POWER SPECTRA 

In this section we show the results on the CMB tem- 
perature and polarization pattern obtained by our mod- 
ifications of the CAMB code. Fig. (JXOJ) shows the var- 
ious contributions to the total CMB temperature and 
polarization angular power spectra from the pure mag- 
netic mode and its correlation with the adiabatic mode. 
Fig. (|TT|) shows the dependence of the total temperature 
power spectrum on the spectral index ris- 
ks is clear from the previous section, the Lorentz force 
of a fully correlated magnetic contribution decreases the 
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FIG. 6: Evolution of baryons velocity for 4 different wavenum- 
bers with (dashed) and without (solid) PMF. k 2 L/pb (dot- 
dashed line) and the solution #j, ate (dotted line) are also plot- 
ted: note how the numerics agree with 6*j, atc at late times. 
The cosmological parameters of the flat AC DAI model are 
Q b h 2 = 0.022, n c h 2 = 0.123, r = 0.04, n s = 1, H = 
72kms~ 1 Mpc" 1 . 
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FIG. 8: Time evolution of baryons (left) and CDM (right) 
density contrasts with vanishing PMF (solid) , fully correlated 
(dashed), fully anti-correlated (dotted) and purely magnetic 
initial conditions (dot-dashed). The other cosmological pa- 
rameters are the same as Fig. JBJl. 





FIG. 9: Time evolution of baryons (dashed) and CDM (solid) 
density contrasts for purely adiabatic with vanishing PMF 
(black) , fully correlated (left panel) and anti-correlated (right 
panel) PMF with vanishing (blue) and non vanishing (red) 
Lorentz force for k — IMpc" 1 . These figures show clearly 
that the Lorentz force and the gravitational contribution are 
of opposite sign, and the Lorentz term is more important. 
The cosmological parameters are the same of Fig. (((J . 




0.0 100.0 1000.0 10000.0 



10.0 100.0 1000.0 



FIG. 7: Evolution of baryons (dotted), CDM (solid) and pho- 
tons (dashed) density contrast for 4 different wavenumbers 
with fully correlated (blue) and without (black) PMF. The 
cosmological parameters are the same of Fig. ([6]). 



density contrasts and therefore the CMB APS in an in- 
termediate range of multipoles before high £ increase. 



VIII. RESULTS FOR THE MATTER POWER 
SPECTRA 



In Fig. (|T3|) we present the results for the linear CDM 
power spectrum evaluated at present time in presence 
of SB of PMF. By analyzing Fourier spectra we have 
checked that the adiabatic results are recovered for k > 



2k d- We compare the results obtained by neglecting or 
by taking into account the Lorentz term. By considering 
the equations evolved and the previous figures, it is clear 
how the Lorentz term treated as in Eq. (f^Bj) is a leading 
contribution for baryons which gives rise to a long-time 
effect as show in Fig. 6. Through gravity CDM is affected 
as shown in Figs. (7-9) and therefore a large feature is 
present in the linear CDM matter PS. 



IX. CONCLUSIONS 

We have investigated the impact of a SB of PMF on 
scalar cosmological perturbations and its impact on CMB 
anisotropies and matter power spectrum. The effects 
on the CMB angular power spectrum is one of the dis- 
tinctive features of stochastic PMF together with non- 
gaussianities and Faraday rotation future missions 
as Planck 1201 will greatly improve the present con- 
straints [iflflq]. 

We have analyzed the SB of PMF in the one-fluid MHD 
approximation Q as a source for cosmological pertur- 
bations and we have inserted such modifications in the 
CAMB cod e Our numerical code improves previ- 

ous studies 1 10L 1 1 111 for the treatment of initial conditions 
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FIG. 10: CMB temperature angular power spectra obtained 
with (B 2 ) — 3 x 10~ 7 Gauss, i%b = — 1, feu = ?r in compar- 
ison with the adiabatic spectrum with vanishing PMF (solid 
line): TT, EE, TE are displayed in the top, middle, bottom 
panel, respectively. The purely magnetic, correlation, fully 
correlated, fully anti-correlated and uncorrelated spectra are 
represented as triple dotted - dashed, dashed, dotted, dot - 
dashed and long dashed lines, respectively. The other cosmo- 
logical parameters are the same as Fig. JB). 
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FIG. 11: In the top panel variation of the CMB temper- 
ature angular power spectrum with tib in comparison to 
the case with vanishing PMF (solid line). In the top figure 
\J (B 2 ) — 8 x 10~ 7 Gauss, ko = 2ty and fully correlated ini- 
tial conditions are considered. The spectral indexes plotted 
are ns = —3/2,-1,1,2 (dotted, dot-dashed, dashed, long- 
dashed lines, respectively). In the bottom panel, variation of 
the CMB angular power spectrum with ko in comparison to 
the case with vanishing PMF (solid line). In the bottom fig- 
ure y/ (B 2 ) = 3 x 10~ 7 Gauss, ub = —1 and ko = 2-7T , ty , 7r/2 
(dotted, dot-dashed, dashed, respectively). In both panels 
the initial conditions are fully correlated and the other cos- 
mological parameters are the same as Fig. (B). 



and exact convolutions for the PMF energy-momentum 
tensor. Note that the present constraints [13, EB] used 
neither the correct initial conditions nor the correct con- 
volutions for the PMF energy density and Lorentz force 
power spectra. Ref. [ll[ uses the correct initial condi- 
tions, but a power spectrum for the PMF energy -density 
with a spectral index which is twice the one for the power 
spectrum of the magnetic field. We have shown exten- 
sively in Sect. VI and Appendix A,B that this is not the 
case. 

We have also shown how the Lorentz term for baryons 
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FIG. 12: In the left panel CMB temperature power spectrum 
obtained with fully correlated PMF with (dashed line) and 
without (dotted line) Lorentz term in comparison with the 
vanishing PMF (solid line). As is clear from the previous 
section, the Lorentz force of a fully correlated magnetic con- 
tribution decreases the density contrasts and therefore there 
is range in which the CMB TT APS is decreased respect to 
the adiabatic case. In the right panel the same figure with 
uncorrelated spectra. In the figures sJ(B 2 ) = 3 x 10~ 7 Gauss, 
ko = 2n and ns = 2 are considered. The other cosmological 
parameters are the same as Fig. (HJ). 
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k (V 1 Mpc"') <B 2 >' /2 =3 10~" n = 2 



FIG. 13: Linear cold dark matter power spectrum obtained 
with fully correlated PMF with (dashed line) and without 
(dotted line) Lorentz term, with uncorrelated PMF and the 
Lorentz force (dot-dashed line) in comparison with the vanish- 
ing PMF (solid line). In the figure y/{B 2 ) = 3 x 10" s Gauss, 
ko = 2-7T and tis = 2 are considered. The other cosmological 
parameters are the same as Fig. ((6]). 

in the one-fluid plasma description Q may lead to a long- 
time effect which we have described analytically in Eq. 
(26). This last point deserves further investigation. 
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X. NOTE ADDED 

While this paper was about to be completed, an arti- 
cle [2l| which computes numerically the convolution in- 
tegrals by taking into account the angular part appeared. 
Ref. [1H does not display the dependence on k of these 
convolution integrals and therefore we cannot compare 
our analytical expressions with their results. 



APPENDIX A: ENERGY DENSITY 

The convolution which gives the magnetic energy den- 
sity spectrum is given in Eq. (|16[) with the parametriza- 
tion for the PS for the magnetic field given in Eq. (|17p . 
Let us compute the energy density of magnetic field ap- 
plying the sharp cut-off to both spectra in the convolu- 
tion according to Eq. (jV A[) . In this case we have two 
conditions of existence to take into account: 

p<k D , \k~p\<k D (Al) 

This second condition poses a fc-dependence on the an- 
gular integration domain and, together with the first 
one, allows the energy power spectrum to be defined for 
< k < 2ko- For simplicity of notation we normalize 
the Fourier wavenumber to k]j and perform the integra- 
tion with this choice. The double integral (over 7 and 
over p) we have to compute must therefore be splitted in 
three parts depending on the 7 and p lower and upper 
/c-depending bounds. This splitting is well displayed in 
(k,p) plane as showed in Fig. (|14j) : in region a the angu- 
lar integration has to be done between —1 and 1, while 
within b and c regions between (k 2 +p 2 — l)/2kp and 1. 
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A sketch of the integration is thus the following: 



1) < fc < 1 

rl-k pi 

dp / drf 

'0 J-l 

2) 1< fc < 2 

f 1 /•! 



dp f 

i-* ^^^^ 



c?7 



i-fc 



dpl a (p,k)+ dpl b {p,k) 



i-fc 



dp 



fe-i 



c?7 



dpl c {p, k) 



fc-i 



(A2) 



The angular integrals can be performed as: 



4 = p 



„n+2 



2- 



fe 2 (l- 7 2 ) 



(fc 2 + p 2 - 2kp 1 ) n/2 d 1 



k 2 + p 2 — 2kpj 

2p n ~ 1 
kn(2 + n)(A + n) 

-|fc - p\ n+2 (k 2 + k(2 + n)p + (1 + An + n 2 )p 2 



(fc + p) n+2 (fc 2 - k{2 + n)p + (1 + An + n 2 )p 2 ^ 



(A3) 



Ib=Ir. = 
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k 2 +p 2 - 
2kp 



P 



p n+2 
.n-1 



2- 



fc 2 (l- 7 2 ) 



(fc 2 +p 2 - 2fcp 7 )" /2 d 7 
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4 + 2n - 8fc 2 n + 6fc 4 n + n 2 



2fc 2 n 2 



fc 4 n 2 

4 i „2 4 
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-16fc 2 p 2 + 24np 2 - \2k 2 np 2 + 6n 2 p 2 - 2k 2 n 2 p 2 + 8p 4 + 6np 4 + nV 
-8|fc - p\ n+2 (fc 2 + fc(2 + n)p + (1 + 4n + n 2 )p 2 



(A4) 



3.2 0.4 . G 



*<l/2 / ( i_fc)*- 



with p > k 
with p < k 
with p > k 



FIG. 14: Integration domains in (fc,p) plane 



Note that the divergent terms at the denominator n 
and n + 2 simply means that the above formulae are not »i 

applicable for n = and n = —2 (logarithmic terms / dp—>(l<k<2 f(k—i) ^P--- with p < fc 

appear in these cases). "'( fe - 1 ) 

A special care must be taken in the radial integral. In 
particular the presence of the term \k — p\ n+2 in both 

integrands, when considering odd spectral indexes, :t is important to study some relevant behaviour of the 
makes necessary to divide the two cases:p < fc and integrands in p. For p ~ 0: 
p > fc.This leads to a further division of the integration 
domain. The scheme of the radial integration is then: 




Jq dp... with p < k 
J*/ 1 fe ' ) dp... with p > k 
Jg 1 ^ <ip... with p < k 



-k n p n+2 , (A5) 
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For p ~ fc the above integrands behave as 

2fc™- 2 



h 



n(n + 2)(n + 4) 
fc™~ 2 



2 n+2 fc n+4 7 



(n + 3) - ((fc -p) 2 ) 2 (n+ l)(n + 4)fc 2 



An(n + 2)(n + 4) 



n(4(n + 4)fc 2 + n + 2) -8((fc-p) 2 ) 2 (re + l)(n + 4)fc 



(A6) 



(A7) 



For I c , p ~ fc cannot be obtained since 1 < fc < 2. It is important to stress that for n > —3 the divergences in p ~ 
and p k are integrable. The coefficients of both leading terms are proportional to fc". 

Following the scheme (IA2[) we can perform the integration over p. Our exact results are given for particular values 

of flfl. 
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APPENDIX B: LORENTZ FORCE 

In order to obtain the complete estimate of the contri- 
bution of PMFs to the perturbation evolution is neces- 
sary to solve the convolution for the Lorentz Force power 



spectrum. The anisotropic stress can be obtained directly 
from its relation with the Lorentz force and the magnetic 
energy density. The convolution which gives the Lorentz 
force is given in Eq. (|23p with the parametrization for 
the PS for the magnetic field given in Eq. (fT?)) . 
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